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We theoretically study Langevin systems with a tilted periodic potential. It has been known that 
the ratio Q of the diffusion constant D to the differential mobility is not equal to the temperature 
of the environment (multiplied by the Boltzmann constant), except in the linear response regime, 
where the fluctuation dissipation theorem holds. In order to elucidate the physical meaning of 
far from equilibrium, we analyze a modulated system with a slowly varying potential. We derive a 
large scale description of the probability density for the modulated system by use of a perturbation 
method. The expressions we obtain show that Q plays the role of the temperature in the large 
scale description of the system and that Q can be determined directly in experiments, without 
measurements of the diffusion constant and the differential mobility. Hence the relation D — fidO 
among independent measurable quantities D, Hd and can be interpreted as an extension of the 
Einstein relation. 
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I. INTRODUCTION 

Technological development in both the manipulation 
and observation of objects on small scales has led to 
a new understanding of the behavior exhibited by me- 
chanical systems of length scales ^ 10^^ m, force scales 
~ 10~^^ N, and the time scales ~ 10"'^ sec. One such 
system that has been studied consists of a small bead 
suspended in a fluid. Because of their simplicity, such 
systems are particularly useful for studying topics rele- 
vant to fundamental physics. Indeed, recently, the fluc- 
tuation theorem |l| and the Jarzynski equality , which 
were derived theoretically as relations universally valid 
for nonequilibrium processes, have been verified experi- 
mentally through experiments on systems consisting of 
small beads 3] and RNA molecules ^] employing optical 
tweezers. 

In studying nonequilibrium systems, we want to dis- 
cover the uniquely nonequilibrium behavior as well as to 
determine what properties of equilibrium systems remain 
even far from equilibrium. We believe that small systems 
are suited for such studies because nonequilibrium effects 
become more significant as the system size decreases. In 
particular, with regard to nonequilibrium systems, we 
are interested in finding new relations between measur- 
able quantities that may be useful in the construction 
of a systematic theory of nonequilibrium statistical me- 
chanics. 

In the present paper, we study the motion of a small 
bead suspended in a fluid of temperature T. The bead is 
confined to move in a single direction, say the x-direction, 
and is subjected to a periodic potential U{x) of period £. 
Such a system can be realized experimentally as a scan- 
ning optical trap system 5], for example. Further, a flow 
with constant velocity can be used to apply a constant 
driving force / to the bead. In this way, it is possible 
to experimentally realize nonequilibrium steady states 
(NESS) for such a bead system. 



The quantity we investigate is the ratio Q of the diffu- 
sion constant D to the differential mobility /id for the 
bead in NESS. In the linear response regime, the ra- 
tio O is identical to the temperature of the environment 
(multiplied by the Boltzmann constant). This relation is 
equivalent to one form of the fluctuation dissipation the- 
orem (FDT). However, because the FDT does not hold 
for NESS far from equilibrium, O is not identical to the 
temperature of the environment. Nevertheless, in the 
system considered there, we find that Q plays the role 
of the temperature in the description of the large scale 
behavior of the system and that & can be determined 
experimentally in a direct manner, without the need to 
measure D and /id. We obtain this result by employing a 
perturbation method to derive the large scale description 
of the probability density. 

II. MODEL 

We assume that motion of the bead is described by the 
one-dimensional Langevin equation 



7i = -^ + /+V^eW, (1) 

where ^(t) is a Gaussian white noise with zero mean and 
unit dispersion. The Boltzmann constant is set to unity. 
We note that there are many physical examples that are 
described by this equation f|. Here, we make two re- 
marks on the form of QJ. First, inertial effects are con- 
sidered to be negligible, because a typical value of the re- 
laxation time of the particle velocity, which is estimated 
to be ^ 10~^ sec 0], is much shorter than the character- 
istic time scale of the phenomena observed in the type of 
experiments we consider. Second, T in is assumed 
to be the temperature of the environment. Although 
the validity of this assumption is not known for NESS 
in general, our result does not depend on the physical 



interpretation of T, because we do not need to use the 
value of T in our analysis. 

The probability density for the position of the particle 
p{x,t) in this system obeys the Fokker-Planck equation 



dp Id 

dt 7 dx 



dx ' I dx 



The steady state density Ps{x; f) is obtained as 



1 



with 



Ux) 



dye 



'l3U(x)+liU{y+x)-l3fy 



(2) 
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where [3 — 1/T and Z is a normalization factor by which 



dxps{x; /) = e. 
The steady state current Vs{f) is derived as 

«s(/) 



l\f,dxl.{x) 



(5) 



(6) 



The exact expression of the diffusion constant D{f) has 
been derived recently for the system under consideration 
13 . By using a different method (See sections IIVI and 
IVI|) . we derive the following form for D{f): 



T \f,dx{I^{x)fl+{x) 



(7) 



where 



I+{x) = / dye 



^l3U(x)-l3U{x-y)-l3fy 



(8) 



Note that is equivalent to the expression derived in 
Ref. 9], which is obtained by simply exchanging and 
/_ in 0. 

In the inset of Fig. we display an example of D{f) 
for the case U{x) = ?7o sin 27rx/^. It can be seen that 
the diffusion is enhanced around fi/T = 2ttUo/T. This 
effect was first reported in Ref._[lCl|| and was subsequently 
analyzed more quantitatively 



III. QUESTION 

In the linear response regime, the mobility /i, defined 
as fi = hm/^o ^^s(/)//, is equal to D(/ = 0)T. This 
is one form of the FDT. However, for NESS far from 
equilibrium, the FDT does not hold. In fact, far from 
equilibrium, there does not even exist an FDT involving 
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FIG. 1: e/T as a function of f£/T in the case that U{x) = 
C/o sin 27r2:/^ with Uo/T = 3.0. e/T = 1 at / = 0, and 
9/r > 1 for / > 0. The inset displays D'y/T as a function of 
f£/T calculated using 



the differential mobility /id = dvs/df , drived from |(HJ as 
follows: 



l)jQdxI-{x)I+{x) 

7 



{lSodxI^{x) 



(9) 



In order to determine quantitatively the extent to 
which the FDT is violated, we define 8 as the ratio of D 
to /id: 



e(/) . 

Md(/) 



(10) 



As displayed in Fig. ^ although the dimensionless quan- 
tity Q /T is unity in equilibrium, it depends on jljT and 
on the form of the periodic potential V (x) when the sys- 
tem is far from equilibrium. In addition to the fact that 
the FDT does not hold far from equilibrium, no means 
of measuring experimentally without measuring D and 
/id has been proposed. Hence, there is no known method 
of deducing /td from experimentally measured values of 
D. While the fluctuation theorem reported in Ref. 
which is valid for NESS far from equilibrium, can be re- 
garded as an extension of FDT > it seems difficult to 
connect it to a physical interpretation of HlOfl . 

Despite the apparent difficulties described above, we 
attempt to propose a method of measuring experi- 
mentally without measuring D and /id. We first note 
that e/T (^ 1 in NESS) is interpreted as the FDT vio- 
lation factor. Recently, stimulated by a proposal of the 
thermodynamic measurement of the FDT violation factor 
in spin glass systems as the " effective temperature" (l^ , 
the feasibility of physical measurements of the effective 
temperature has been investigated in numerical experi- 
ments modeling a sheared glassy material jl3| . a dri ving 
system near jamming jl4| , and driven vortex lattices jlSj . 
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Among them, Berthier and Barrat have proposed a novel 
measurement method of the FDT violation factor as the 
effective temperature in shared glassy systems [13l ■ They 
have used a tracer particle of large mass as an effective 
thermometer and demonstrated that the kinetic energy 
of the tracer is related to the FDT violation factor. This 
result exhibits a clear relation between the FDT violation 
factor and the measurable effective temperature. 

Here too we seek to determine whether there is a realm 
in which Q can be interpreted as the effective tempera- 
ture in NESS and, hence, whether it can be measured 
independently of D and /id- 



IV. RESULT 

We propose a method of measuring Q by adding a 
slowly varying potential which plays the role of an ef- 
fective thermometer to the system under consideration. 
In order to demonstrate this method, we study with 
U{x) replaced by a potential U{x) + V{x), where V{x) is 
a slowly varying periodic potential V{x) of period L ^ £. 
We investigate the behavior of the probability density on 
length scales larger than L. In order to make the sepa- 
ration of scales explicit, we define e = £/L and the large 
scaled coordinate X = ex. The probability density in 
this system obeys the Fokker-Planck equation 



dp 
'di 



ld_ 

7 dx 



— +e— Ap + T— 
dx dX j dx 



(11) 



where we have defined V{X) = V{x). 

We extract the large scale behavior of p{x, t) by intro- 
ducing a slowly varying field Q{X,t) as 



p{x,t) = Ps{x; frn) [ Q + ea{x;fm)^ + e^b{x;f^)^ ^ 



+ ec{x;U)^Q + Ois')], 



dX^ 
(12) 



where fm = f — edV/dX. Because /m is a function 
of X = £x, the X dependence of Ps{x; /,n) appears in the 
two ways, as an explicit dependence and as a dependence 
through /m. Note that Ps{x; /m) is a periodic function in 
the sense that ps{x + fm) = Ps{x;fm)- The functions 
a(a;;/m), b{x] fm) and c{x] fm) are similar; that is, they 
are periodic functions of x in the same sense and depend 
on /in. Their functional forms are determined below. 

Substituting ifT^ into ((TT|l . we obtain 



dQ_ 

dt 



eA 



dQ 
dX 



(13) 



where A, B and C can be expressed by 
A = 



U-^pAp-^{l + a') + -a\ (14) 
V 7 / 7 



B 



T 
—1 

7 



-Ps P 



-{1 + 2a' + b") ~ aA, 

7 



T 



Ps'C 



T 



7 O/m 



(15) 



(16) 



where ' represents the partial derivative with respect to 
X, that is, a'(x]fm) = da(x] f^)/dx and so on. Then 
we can choose a, b and c so that A, B and C do not 
depend on x explicitly but depend on x through the X 
dependence of /m. After a straightforward calculation, 
we find that, subject to this condition, A^ B and C are 
uniquely determined as 



-4 = -Ws(/m), 

B = DiU, 



C 



dfn 



(17) 
(18) 

(19) 



(The derivation will be presented in section IvTI ) Using 
this result, we rewrite as 



dQ 
dt 



d 

'dX 



-Vs{fm)Q + D{U 



dQ 

dx 



(20) 



Recalling that /m — f ~ edV/dX, we see that this equa- 
tion is equivalent to the Fokker-Planck equation 



dQ I d 
'dt ~ f^dX 



dX 

where we have defined 



e^-F]Q + ee?^+Ois^] 



dX 



r(/) EE A*d(/)-\ 

F{f) EE V,{f)flMy\ 



(21) 

(22) 
(23) 



For reference, in Fig. |21 we present graphs of F and F as 
a function of / for the model considered in Fig. ^ 

We note that H21|l has the same form as with 
the parameters F, F and O corresponding to 7, / and 
T, respectively. Thus, replacing {'j,f,T,U{x)) in JHl by 
iT,F,Q,V{X)), we obtain the current for NESS of this 
modulated system. From this result, by measuring the 
steady state current for several forms of V{x), we can de- 
termine the values of F, and experimentally, where 
we note that these values do not depend on the choice of 
V{x). This fact implies the existence of the new relation 
D ~ fi(i(^ among the independent measurable quantities 
D, and O in NESS far from equilibrium. This new re- 
lation D = ii^Q can be thought as the extended Einstein 
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FIG. 2: The left axis represents r/7, which is displayed as 
a function of /^/T (solid curve), and the right axis rep- 
resents Fl/T, also displayed as a function of fl/T (dash- 
dotted curve), in the case that U(x) = Uosm2-KX / with 
Uo/T = 3.0. The dotted line corresponds to F = f . It is 
seen that in the linear response regime, the effective force F 
is equal to /. 



relation. Furthermore, from the correspondence between 
O in H21() and T in it is evident that Q plays the role 
of the temperature for the large scale behavior of the sys- 
tem. In this way, we have arrived at the main claim of 
this paper, the physical interpretation of Q. 

Here we address two remarks on the main claim. First, 
one may naively expect that on the large scale can be 
effectively described by 



x = v,{f) + ^2D{r)(,{t). 



(24) 



Although such an effective description is valid, we em- 
phasize that the effective temperature is not determined 
from this description. In order to have the correspon- 
dence with (|21|l . we need one more quantity in addition 
to D and Ug. In our analysis, by adding the slowly vary- 
ing potential to the system, F, F and O are determined. 
We thus interpret the slowly varying potential as an ef- 
fective thermometer for the original system described by 

Second, we have another method of measuring Q by 
considering the case in which a slowly varying potential 
V{e{x — Vst)) that moves with constant velocity Vg re- 
places V{ex) in Hll|) . Although the resulting system is 
more complicated than that considered above, the large 
scale behavior in this case is actually simpler. Indeed, 
using the same method as above, we obtain the following 
equation for Q{Y,t) describing the large scale behavior 
of the system: 



dQ(Y,t) _ I _d_ dVjY) 
di ^ f^W \ ^ dY ' 



(25) 

Here, we have introduced the large scaled moving coor- 
dinate Y = e{x — Vgt). This equation is identical to the 




FIG. 3: (V^) /T versus A/T obtained in the numerical ex- 
periment. The dotted line represents the fitting curve of the 
form ^ in which e/r= 1.67(4). 



Fokker- Planck equation describing the time evolution of 
the probability density in an equilibrium state with tem- 
perature G. Therefore, for example, O is obtained by 
measuring the statistical average of V{Y). 



V. NUMERICAL EXPERIMENT 

Let us demonstrate that the measurement method of 
Q by use of a slowly varying potential works well in 
numerical experiments. We study numerically JQ) with 
periodic boundary conditions for the case that U (x) = 
t/p sin 27ra;/£ and that the system size is L. The param- 
eter values are chosen as follows: L/£ = 50, Ua/T = 3 
and fi/T = 16. 

We measure the effective temperature for this system 
by using a moving potential V{y) in the form 

V{y) = 2Ay/L for 0<y<L/2 

= 2A-2Ay/L for L/2<y<L, (26) 

with y = mod(x — Vgt, L). The statistical averages of V 
in NESS, {V), are measured for several values of A. As 
displayed in Fig. |21 the graph (A, {V)) is fitted well to 
the form 



A 



exp(A/e) - 1' 



(27) 



which can be expected from (|25|l . Using this fitting for 
the experimental result, we evaluate the value of Q/T 
as 0/r = 1.67(4). This value should be compared with 
the theoretical one, 1.677, which is calculated from H1(J|I . 
Thus we conclude that our measurement method of Q 
works well in experiments. 



5 



VI. TECHNICAL DETAILS 

We regard (|14() as an ordinary differential equation 
(ODE) for a{x; /,„) under tlie condition that A does not 
depend on x. This ODE has periodic solutions of x 
only when A satisfies a certain condition, which provides 
uniquely A. 

Let us solve the ODE We first define 

d{x; fni) = a'{x; /m)ps(x; /m)- (28) 
Then 1)14(1 becomes 

T T 

—d' - v^a' -Vs + -pi - Aps = 0. (29) 

7 7 

Integrating this equation over the range [0, £], we obtain 

A^-v,, (30) 

where we have used a requirement that Td/^ — Vga is a 
periodic function. 

Under the condition 1(30(1 . we can derive periodic solu- 
tions a(a;;/,n). The integration of 129(1 leads 

T T 

—a'ps — Vga — VgX H Ps + VgH + VgKi — 0, (31) 

7 7 



where Ki is a constant and we have defined 



H{x]fn^) = / dypsiyjr^). 
Jo 



(32) 



Substituting the expression 

a(x; /m) = H{x; /,„) - a; + a(x; /m) (33) 
into ((^ . we rewrite l(M|l as 



T T 

— Psa' - Vsd H _Ps + VsKi = 0. 

7 7 



(34) 



e---' — I , (35) 



Noting the relation 

Vs (pse 
we solve ((34(1 as 

d{x;U) = -Ps{x;fmr'e-P^'"^^^~f-^^^x;U 

+ i^2ft(x;/„,)-ie-'5(^(^)-^»^)+i^i, (36) 

where we have defined 

H^; /m) = r dyp^iy; f^^^)^e^myJn.)~Uy) _ (37) 



From the condition a(0; /m) = a(€; /m), the constant 
is determined as 



1 



1 _ e-/3/m^ 



$(^;/m). 



(38) 



((33(1 . 1(36(1 . 1(3 7() and 1(38(1 provide all periodic solutions of 
ODE ((Tljl . 

Next we study (115(1 by repeating the similar analysis. 
Defining 

Hx; /m) = b'{x; fm)ps{x; /m), (39) 
we rewrite ((15(1 as 
T T T 

—h' — Vgb' — Vga H (1 + 2a')ps H ap^ — Bps + aVsPs = 0. 

7 7 7 

(40) 

Integrating this equation over the range [0, £], we obtain 

B=- + — / dxa'p, - ^ / dxa{l - p,), (41) 
1 l*- Jo * Jo 

where we have used a requirement that Th/^ — Vsb is 
a periodic function of x. We can derive the expression 
of B by substituting a{x; fm) into (|^ . However, this 
expression is very complicated. We now simplify this. 

Eliminating a' in (|IT|l by use of (|^ and ((M)l . wc ob- 
tain 

B^^ dx{x-H) + ^ dxpsa-VsKi. (42) 
* Jo ^ Jo 



We can simplify this equation as 



B = dxpstt — VsKi 



where we have used the equality 



I dx{x - H){1 - ps) = 0. 
Jo 



(43) 



(44) 



Substituting (|!^ into (|^. we obtain 



B 



£{1 - e-/3/»^) Jo 
m£;U-{l-e-^f-'Mx;U]. (45) 



Here we have an identity for an arbitrary periodic func- 
tion (l){x) with the period £: 



dy(j){y)e-l^f"^y -{\-e-^f-^^) / dy(j){y)e-^f-^y 
lo Jo 

^ Q-fiUx dy(j){y + x)e-'^f^^y . (46) 
Jo 

Putting (j) = p^e^^ in ((46(1 . we simplify ((45(1 as 



B = 



(1 _ e-/3/»0^ 



X / d2/(Ps(2/ + x;/„,))2e^^(^+^)-''/»^. (47) 
Jo 
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Using Q and 10), we can derive 
B = D{U) 

with 

Finally, defining 



(48) 



g{x; /m) = c'{x; fm)Ps{x; /m), (49) 
we rewrite (fT^ as 

T , I dvs T dp'g 



-g -VsC . 
7 a/m 7 oiu 



PsC = 0. (50) 



Integrating this equation over the region [0, £], we obtain 



(51) 



where we have used a requirement that Tg/7 — UgC is a 
periodic function of a;. 



VII. DISCUSSION 

In conclusion, we have proposed the method of measur- 
ing the effective temperature of the Langevin system ^ 
by using a slowly varying potential and have found that 
this effective temperature is equal to defined by (|10|l . 
The independence of measurements of the quantities 

and O makes us interpret D — fidQ as the extended 
Einstein relation of the Langevin equation. This signifi- 
cant result was obtained by analyzing the Fokker-Planck 
equation with a slowly varying potential. 

At the end of this paper, we shall present remarks on 
(|25|l which is derived by the analysis of the system with 
the moving potential. H25() provides two important in- 
sights in addition to the direct measurement method of 

e. 

The first insight obtained from H25() is related to the 
interpretation of the extended Einstein relation D — ^d© 



among independent measurable quantities. Because H25|) 
is identical in form to equations that describe equilibrium 
systems, D = Q/T (the Einstein relation in the linear 
response theory) should hold. From 122|) . this Einstein 
relation yields D — /j^Q. As is well known, the Einstein 
relation is closely connected to the existence of detailed 
balance for fluctuations. However, note that fluctuations 
described by with / 7^ do not satisfy the detailed 
balance condition, as can be easily checked by using the 
steady state density ©. Therefore, we find that the de- 
tailed balance condition is recovered through the coarse- 
graining procedure yielding and that ((Tn|l can be 
understood as the result of the recovery of detailed bal- 
ance with respect to the canonical distribution for the 
temperature 8. 

The second insight obtained from H25|l is related to the 
extension of thermodynamics to NESS. When we assume 
that V(Y) takes a tanh-like form with amplitude A (that 
is, V{oo) — T^(— 00) = A), we have 



Q+-Q- 




(52) 



with Q± = \iuiY^±(x> Q{Y)- Then, the chemical poten- 
tial extended to NESS can be defined usinK_p3 in a sim- 
ilar way as the case of driven lattice gas [l^ ■ Therefore, 
it may be possible to incorporate the idea of the effective 
temperature into a theoretical framework of thermody- 
namic functions extended to NESS. A study with this 
aim treating a wide class of nonequilibrium systems, in- 
cluding many body-systems, is now in progress. 
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